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We propose a new method for calculating optical defect levels and thermodynamic charge- 
transition levels of point defects in semiconductors, which includes quasiparticle corrections to the 
Kohn-Sham eigenvalues of density-functional theory. Its applicability is demonstrated for anion 
vacancies at the (110) surfaces of III-V semiconductors. We find the (+/0) charge-transition level 
to be 0.49 eV above the surface valence-band maximum for GaAs(llO) and 0.82 eV for InP(llO). 
The results show a clear improvement over the local-density approximation and agree closely with 
an experimental analysis. 

PACS numbers: 71.15.Qe, 71.45.Gm, 73.20.Hb 



The electrical and optical properties of semiconductors 
depend sensitively on the electronic structure in the gap 
region and can hence be modified dramatically by the 
presence of native defects and impurities that introduce 
unwanted additional states inside the fundamental band 
gap. Most importantly, such electrically active defects 
can trap charge carriers (electrons or holes), counteract- 
ing the effect of intentional doping. If their concentration 
is sufficiently high, this process can lead to a full com- 
pensation of implanted acceptors and donors and thus 
eventually to Fermi- level pinning. Besides, electron-hole 
recombination at deep defects drastically reduces the life- 
time of minority carriers, and transitions involving defect 
states inside the band gap may dominate optical absorp- 
tion. This is especially relevant at surfaces and interfaces, 
where the crystal termination and the contact with other 
phases naturally give rise to a high number of structural 
defects that have been linked to the formation of Schott- 
ky barriers [l[ . The central quantities are the optical de- 
fect levels inside the band gap and the charge-transition 
levels. The former can, in principle, be probed by direct 
(filled states) or inverse (empty states) photoemission. 
The Franck-Condon principle is well justified, as the re- 
arrangement of the atoms happens on a much slower time 
scale than the electron emission or absorption, but the 
coupling to the atomic lattice may be visible in the line 
widths and shapes. The optical defect levels contain the 
full electronic relaxation in response to the created hole 
or the injected electron, however. The charge-transition 
levels, on the other hand, are thermodynamic quantities 
and specify the values of the Fermi energy where the 
charge state of the defect changes. Therefore, they are 
affected noticeably by the atomic relaxation taking place 
upon the addition or removal of an electron. 

Despite considerable efforts, a reliable determination of 
the optical defect levels and the related charge-transition 
levels of deep defects still poses a very difficult challenge. 
Experimental investigations are thwarted by the fact that 



many traditional spectroscopic techniques arc not appli- 
cable due to the low concentration of native point defects, 
while capacitance methods like deep-level transient spec- 
troscopy are very sensitive but provide no elemental 
or structural information to identify the type of defect. 
For surfaces, at least, Ebert et al. [3j now demonstrated 
how the electronic structure of individual defects can be 
deduced using a combination of scanning tunneling mi- 
croscopy (STM) and photoelectron spectroscopy, provid- 
ing the first reliable experimental analysis of a charge- 
transition level for the P vacancy at InP(llO). Theoreti- 
cal approaches, on the other hand, must include accurate 
exchange-correlation contributions, the coupling between 
electronic and lattice degrees of freedom, and, in general, 
require the treatment of open systems in which the num- 
ber of particles is not constant Q . Previous studies that 
employed density-functional theory (DFT) in the local- 
density approximation (LDA) mf indeed suffered from 
fundamental limitations [3, @, 0713] ■ For example, Ebert 
et al. [3] noted "that the systematic error for the calcu- 
lated energies of the charge transfer levels is too large to 
identify the symmetry of the vacancy on the position of 
the defect level only." In order to overcome this problem 
we here propose a new computational approach, broadly 
applicable to defects in the bulk as well as at surfaces, 
that combines DFT with many-body perturbation the- 
ory. As an example, we examine the optical defect levels 
and the thermodynamic charge-transition levels of anion 
vacancies at GaAs(llO) and InP(llO). The results are in 
close agreement with the experimental analysis Q. 

The geometry of anion vacancies at the (110) surfaces 
of III-V semiconductors is now well understood thanks 
to a combination of experimental and theoretical studies. 
For p-doped materials STM images of the filled states un- 
der negative bias feature a localized hole at the position 
of the missing anion surrounded by a voltage-dependent 
depression [9|. The latter is due to a downward local 
band bending, indicating a positive charge of the va- 
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cancy. The charge state is, in fact, established as +1 
[lot , which is also predicted by electronic-structure cal- 
culations 0,0)0)1^]. STM images acquired under positive 
bias probe the empty p z -like orbitals of the cation sublat- 
tice and show an enhancement of the cations surround- 
ing the vacancy, initially wrongly interpreted as an up- 
ward relaxation of those atoms [9( but now understood as 
arising from the local depression of the electron density. 
DFT-LDA calculations actually show that an inward re- 
laxation of the two Ga atoms enclosing the As vacancy 
at GaAs(llO) is consistent with the observations 
[a III- The symmetry of the positively charged anion 
vacancies was initially a matter of controversy [y, [fl] ■ In 
STM images they appear symmetric, but in a combined 
experimental and theoretical study of Vp at InP(llO) 
Ebert et al. [|| explained the observed features as result- 
ing from the thermal flip motion between two degenerate 
asymmetric configurations. This interpretation later re- 
ceived further confirmation [111 ]. In n-doped materials, 
on the other hand, the vacancy is in a charge state of 
— 1 [13], and DFT-LDA calculations predict a symmet- 
ric relaxation for this configuration as well as the neutral 
vacancy 0,0. 

The anion vacancies at GaAs(llO) or InP(llO) give 
rise to three nondegenerate electronic states, labeled la', 
la", and 2a'. While the la' state is located several eV 
below the valence-band maximum and always filled with 
two electrons, and the 2a' state is too high in energy 
to become populated, the la" state lies inside the band 
gap. It is this state, therefore, that is relevant for the 
discussion of charge-transition levels. Depending on the 
level of doping it may be occupied by zero, one, or two 
electrons, which corresponds to the positive, neutral, and 
negative charge state, respectively. 

The formation energy of a surface vacancy with charge 
state g, relative to that of the neutral defect, is given by 
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forn 



(q/0) = E vac (q, Q q ) - £ vac (0, Q ) + qe F , (1) 



where E vac (q,Q) denotes the total energy of a surface 
featuring a single vacancy with the actual electron pop- 
ulation q G {+, 0, — } and geometry optimized for charge 
state Q € {Q + ,Qo,Q-}. The final term accounts for 
the transfer of the charge q between the defect level and 
the electron reservoir, i.e., the Fermi energy ep. The 
charge-transition levels t q l q are defined as the values 
of ep where the charge state of the vacancy changes, 
i.e., where £ form (q/0) = £ ,form (g'/0), and convention- 
ally given relative to the surface valence-band maximum. 
For the systems considered here the interesting transi- 
tions are e + /° and e°/ _ . For example, the former is 



H +/o =£ vac (0,Q o )-£ vac (+,Q + ). 



(2) 



All previous DFT-LDA calculations for surfacepoint de- 
fects evaluated this energy difference directly 0, @, 0, [1] , 
but this approach leads to systematic errors that arise 



because the total energies E™ c (0,Q o ) and E vac (+,Q + ) 
refer to systems with different electron numbers. As is 
well known, the exact exchange-correlation potential in 
DFT exhibits a discontinuity upon addition or removal 
of an electron [13], which is not contained in the LDA 
or other jellium-based functionals. Besides, aspects like 
the self-interaction are treated inappropriately. As a con- 
sequence, the band gaps of semiconductors and the en- 
ergies of localized defect states are not given correctly: 
for the P vacancy at InP(llO) the experimentally deter- 
mined e + /° level of 0.75±0.1 eV Q contrasts noticeably 
with the calculated values 0.52 eV @| and 0.388 eV @|. 
The variation between the two theoretical results can be 
traced to differences in the pseudopotentials and parallels 
the variation of the corresponding band gaps. 

To arrive at a more accurate quantitative method that 
corrects the above-mentioned severe shortcomings of the 
LDA we rewrite Eq. ([2]) as 



[£ vac (+,Qo)-£ vac (+,Q+)] 

+ [£ vac (0,Q o )-£ vac (+,Qo)] 



(3) 



by adding and subtracting the total energy E vac (+, Qq) 
of a system with the geometry of the relaxed neutral va- 
cancy but a charge state q = +1. In this way the charge- 
transition level is decomposed into two separate contribu- 
tions. The first describes the structural relaxation energy 
for the positive charge state and two different geometries: 
that of the neutral and that of the positive charge state. 
It is always positive. As the electron number remains 
constant, the problem of the discontinuity does not arise, 
and DFT-LDA is perfectly applicable. The second term 
equals the ionization energy of the neutral defect, where 
the removed electron is transfered to the reservoir, i.e., 
the Fermi energy. Determining the ionization potential 
from the quasiparticle band structure requires a correc- 
tion of the Kohn-Sham eigenvalues, for which we employ 
many-body perturbation theory. Specifically, we use the 
Go Wo approximation for the electronic self-energy 14 1. 
This approach is known to yield reliable band gaps for 
III-V semiconductors [Hj] and their surfaces [lij ]. 
In the same spirit, e°/~ can be written as 
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[£ vac (-,Q-)-£ vac (- Qo)] 
+ [£ vac (-,Q )-£ vac (0,Qo)] 



(4) 



The first term in this case describes the energy differ- 
ence of the vacancy with q = — 1 between its own equi- 
librium geometry and that of the neutral charge state. 
It is always negative. The second term equals the elec- 
tron affinity of the neutral charge state. In principle, 
the self-energy of the neutral charge state yields both 
the ionization potential and the electron affinity, which 
correspond to the energy of the highest occupied and the 
lowest unoccupied quasiparticle state, respectively. From 
a computational point of view, however, this procedure is 
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inconvenient, because the neutral defect has an odd num- 
ber of electrons and requires a spin-polarized calculation. 
Instead, we follow an equivalent approach and extract 
the energy levels from two separate calculations for non- 
spin-polarized systems with an even number of electrons. 
In practice, we thus determine E vac (0, Q ) - E vac (+, Q ) 
as the electron affinity of the positive charge state and 
E vac (-,Q ) - E vac (0,Q o ) as the ionization potential of 
the negative charge state, both in the Qo geometry. 

In the following we apply the expressions derived above 
to anion vacancies at the (110) surfaces of GaAs and InP. 
To determine the defect geometries we use DFT together 
with norm-conserving pseudopotentials 17] and the LDA 
exchange-correlation functional The surfaces are 

simulated using a supercell with a (2x4) periodicity in 
the [001] and [110] directions, consisting of six atomic 
layers separated by a vacuum buffer equivalent to four 
layers. A single vacancy is created at one side of the slab, 
while the dangling bonds at the other are passivated by 
pscudoatoms with noninteger nuclear charges of 0.75 and 
1.25 for anion and cation termination, respectively. This 
mimics the continuation of the substrate by a III-V bulk 
layer. We use the theoretical lattice constants 5.55 A for 
GaAs and 5.81 A for InP to prevent errors resulting from 
a nonequilibrium unit-cell volume during the surface re- 
laxation. The integration in reciprocal space is carried 
out with a mesh corresponding to eight k-points in the 
two-dimensional Brillouin zone of the (lxl) unit cell of 
the defect-free surface. In the case of charged defects we 
apply a uniform compensating background in order to en- 
sure overall charge neutrality. For GaAs we thus obtain 
the relaxation energies E vac (+, Q Q )-E vac (+, Q+) = 0.30 
eV and E vac {~, Q_) -£ vac (-, Q ) = -0.13 eV. The cor- 
responding values for InP are 0.20 eV and —0.17 eV. 

Before presenting our quasiparticle results, we first 
calculate the charge-transition levels strictly within the 
LDA by invoking the Slater-Janak transition-state ap- 
proach [191 ]. where the ionization potential and the elec- 
tron affinity equal the eigenvalue of the la" level de- 
termined self-consistently with the noninteger occupancy 
0.5 and 1.5, respectively. The transition state corrects, at 
least partially, the erroneous self-interaction of the LDA 
but not the discontinuity problem. The resulting energy 
contributions are displayed in Fig. [T] relative to the sur- 
face valence-band maximum. The occurrence of slightly 
negative values in some cases implies that the defect level 
actually falls below the valence-band maximum. This 
is an artefact of the constrained nonequilibrium geom- 
etry: if the atomic structure is allowed to relax, then 
the defect level always lies inside the band gap. For all 
systems studied here we find that the transition-state ap- 
proach yields the same results, with a deviation of less 
than 0.01 eV, as a straightforward evaluation of Eq. ([2|) 
and the corresponding formula for e°/~ within the LDA. 
Quantitatively, our calculated value of 0.47 eV for e + /° in 
InP exhibits the same systematic underestimation of the 
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FIG. 1: Position of the optical la" defect level of the anion 
vacancies at GaAs(llO) and InP (110) as a function of the oc- 
cupation number. The geometry is identical in all calculations 
for the same material and optimized for the neutral charge 
state. The Slater-Janak transition-state approach (filled cir- 
cles) yields the ionization potentials and electron affinities in 
the LDA. More accurate results are obtained by calculating 
the quasiparticle corrections within the Go Wo approximation 
(filled squares). In either case the energy zero is set to the 
respective surface valence-band maximum. 



experimentally derived charge-transition level as earlier 
studies at this level of approximation 0, Q . 

In order to determine the electronic contribution to 
the charge-transition levels more accurately we employ 
many-body perturbation theory. The energies derived 
within this framework correspond directly to the values 
measured in direct or inverse photoemission. We follow 
the usual approach to calculate the quasiparticle energies 
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+ (<P&Wla»)-V xc \ V ®,) 



(5) 



as a first-order correction of the Kohn-Sham eigenvalues 
e^». Here £ is the complex, nonlocal, and frequency- 
dependent self-energy, which we evaluate in the GoWo 
approximation using the Green function Go of the un- 
derlying Kohn-Sham system. Our numerical implemen- 
tation is based on the space-time method [20| . The lo- 
cal exchange-correlation potential V xc must be subtracted 
from the self-energy to avoid double counting. A more 
detailed discussion of our computational method can be 
found in Ref. 2l|. The quasiparticle corrections, but not 



the Kohn-Sham eigenvalues, are obtained from a smaller 
(2x2) surface cell, which reduces the computational ef- 
fort considerably. Although we investigate charged sys- 
tems, the relative self-energy shifts are insensitive to the 
size of the supercell, because they only include exchange- 
correlation effects and no electrostatic Hartree contribu- 
tion. Keeping the atomic positions fixed at the opti- 
mized geometry for the neutral vacancy, we performed 
separate GqWq calculations for the positive and negative 
charge states. The calculated la" single-particle energies 
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TABLE I: The charge-transition levels associated with the As 
vacancy at GaAs(llO), given in eV. Values in brackets refer to 
the constrained symmetric relaxation of the positively charged 
vacancy. The quasiparticle band gap of 1.55 eV in this work, 
calculated at the theoretical lattice constant 5.55 A, is close 
to the experimental value 1.52 eV. 
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e o/- 


LDA (this work) 


0.24 


(0.07) 


0.15 


LDA (Ref. [6]) 


0.32 




0.4 


LDA (Ref. [7j) 




(0.10) 


0.24 


Go Wo (this work) 


0.49 


(0.32) 


0.60 



TABLE II: The charge-transition levels associated with the P 
vacancy at InP(llO), given in eV. The quasiparticle band gap 
of 1.52 eV in this work, calculated at the theoretical lattice 
constant 5.81 A, is close to the experimental value 1.42 eV. 



LDA (this work) 0.47 0.54 

LDA (Ref. [3]) 0.52 

LDA (Ref. [§]) 0.388 0.576 

G W (this work) 0.82 1.09 



Expt. (Ref. [3]) 0.75±0.1 



£ia"(+, Qo) an d £ia"(~ , Qo) ar e shown in Fig.[T]with and 
without the self-energy correction. 

In Table U we summarize the results for the As vacancy 
at GaAs(llO). Values in brackets refer to the constrained 
symmetric relaxation of V^, and are included for the pur- 
pose of comparison with earlier studies. In contrast to 
Refs. 0, 0], which found a stable neutral charge state 
within a narrow energy window, our own calculation at 
the level of the LDA indicates e + /° > e°/ _ and hence a 
direct transition from the positive to the negative charge 
state, but the small energetic separation is within the 
uncertainty of the calculation. With the Go Wo approx- 
imation we find a reversed ordering, which implies the 
existence of a stable neutral charge state, and a slightly 
increased splitting of the charge-transition levels. 

The charge-transition levels for the P vacancy at 
InP(llO) are listed in Table El Our LDA results are 
similar to those reported previously 0, @| and well below 
the experimentally deduced value of 0.75±0.1 eV Q . The 
Go Wo approximation, on the other hand, yields a value 
for e + /° that lies within the experimental error bar. 

In conclusion, we have developed a general computa- 
tional scheme for the optical defect levels and thermody- 
namic charge-transition levels of point defects in semicon- 
ductors. The method is broadly applicable to the bulk 
as well as to surfaces. It relies on a separation of struc- 
tural and electronic energy contributions that can be ac- 
curately evaluated within DFT and many-body pertur- 
bation theory, respectively. In this way the discontinu- 
ity of the exchange-correlation potential as well as other 



shortcomings of the LDA are treated appropriately. Our 
calculated (+/0) charge-transition level for the P vacancy 
at InP(llO) is in close agreement with the experimental 
analysis, confirming the accuracy of this method. 
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